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Abstract 

Voyager 2 observations revealed that Neptune’s internal luminosity is an or¬ 
der of magnitude higher than that of Uranus. If the two planets have similar 
interior structures and cooling histories, Neptune’s luminosity can only be 
explained by invoking some energy source beyond gravitational contraction. 
This paper investigates whether Centaur impacts could provide the energy 
necessary to produce Neptune’s luminosity. The major Endings are (1) that 
impacts on both Uranus and Neptune are too infrequent to provide lumi¬ 
nosities of order Neptune’s observed value, even for optimistic impact-rate 
estimates, and (2) that Uranus and Neptune rarely have significantly differ¬ 
ent impact-generated luminosities at any given time. Uranus and Neptune 
most likely have structural differences that force them to cool and contract 
at different rates. 
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1. Introduction 


While the ice giants may have similar interior structures (e.g. Podolak et 


1995; Fortney and Nettelmann, 2010), their internal luminosities differ 


by a factor of 10. From Voyager 2 IRIS radiometer observations. Pearl and 


Conrath (1991) calculated an internal luminosity of log L/Lq = —11.024 for 


Neptune, while Pearl et ah (1990) found an internal luminosity of log L/Lq = 
— 12.054 for Uranus. The 2.7M^ mass difference between the two planets is 
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not enough to explain the luminosity difference: the internal power generated 
per unit mass is 3.22 x 10“^ erg g“^ s“^ for Neptune and 3.92 x 10“^ erg g“^ s“^ 
for Uranus (Pearl et ah, 1990; Pearl and Conrath, 1991|). Multiple theories 


explaining the energy balance of the ice giants have been put forward, in¬ 


cluding stable stratification in Uranus’ interior (Podolak et ah, 1990), early 


and efficient heat transport by baroclinic instability in Uranus (Holme and 


Ingersoll 

1994), 

and efficient capture of strongly interacting dark matter by 

Neptune 

|Mitra 

2004; 

Adler 

2009 

)• 


One energy source that has not been investigated in connection with ice 
giant energy balance is impact heating. Given a sufficient supply of Cen¬ 
taur^ impacts onto Neptune could be frequent enough to boost Neptune’s 
luminosity to observed values. Indeed, meteoroid impacts onto the moon 
generate hashes of optical light, first observed by Dunham et ah (1999). 
Energy deposited in ice giant atmospheres by Centaurs that penetrate the 
photosphere would not be released instantly to space, as in the case of lu¬ 
nar meteoroid hashes, but would instead be radiated away on a ~ 100-year 
timescale (Conrath et ah, 1990). This paper explores the possibility that 


Centaur impacts may contribute significantly to ice giant luminosities. 

The investigation begins with an order-of-magnitude calculation of the 
typical Centaur size required to produce Neptune’s luminosity with impacts 
alone, treating impacts as a steady-state process. Next, we explore different 
impact rates and break the steady-state assumption, treating impacts as 
a stochastic process. Section contains estimates of the total number of 
Centaurs, which we use as a scaling factor for published impact rates. Section 
[^describes a Monte Carlo approach to computing a cumulative probability 
distribution of planet luminosity. Results and conclusions are presented in 
Section |5l 


^While we use the word “Centaur” loosely to describe any object that may impact 
Uranus or Neptune, Jewitt ( |2009 ) defines Centaurs as comets whose dynamics are con¬ 
trolled by perihelion and/or aphelion interactions with giant planets, such that perihelia 
q and semimajor axes a are in the range 5.2 < (g, a) < 30.0 AU. The Minor Planet Center 
website, minorplanetcenter.net/blog/asteroid-classification-i-dynamics/, defines a Centaur 
as an asteroid with q > 5.2 AU and a < 30.0 AU. 
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2. Impact-Induced Luminosity: Order-of-Magnitude Estimate 


To get a basic idea of how much impacts contribute to ice giant luminosi¬ 
ties, we assign a constant value M to each planet’s accretion rate and assume 
a constant accretion-generated luminosity. The impact-generated luminosity 
is then 

( 1 ) 


R 


with 


where M is the planet mass and R is the planet radius. Equating Limp 
Neptune’s present luminosity requires M = 4 x 10^”^ g yr“^. Based on simu¬ 
lations of diffusion from the Kuiper Belt to the inner Solar System, [Levison 
and Duncan (1997) found that comets impact Uranus and Neptune slightly 


more than once per thousand years. To deliver the average M quoted above, 
most of the impacting Centaurs with a density of ~ 1 g cm“^ would have to 
have radii over 40 km. Such a large average Centaur size can be ruled out by 


crater observations; for example. Stern and McKinnon (2000) calculated that 


the largest craters detected on Triton were created by impactors with radii of 
1-5.5 km. The occultation surveys of Roques et ah (2006) and Schlichting et 


al. (2009, 2012) also indicate a Centaur/Kuiper Belt Object size distribution 
heavily biased toward sub-kilometer bodies. 


In the steady-state scenario where Ljmp is constant and the Levison and 


Duncan (1997) impact rate applies, impacts clearly cannot drive Neptune’s 
internal heating. Explaining Neptune’s luminosity with impacts alone re¬ 
quires one of two scenarios: (1) a substantially higher impact rate, which is 
possible if [Levison and Duncan (1997) underestimated the total number of 
Centaurs; or (2) a recent giant impact that has driven Neptune’s luminosity 
to an above-equilibrium value. The rest of this paper examines scenarios (1) 
and (2). 


3. Total number and size distribution of Centaurs 

Determining the frequency and energy of impacts on ice giants requires 
knowing both the total number of Centaurs and their size distribution. The 
number of Centaur detections is too small to reconstruct a size distribu¬ 
tion based on observations alone: only 7 Centaurs met the “secure orbit” 
standards used by the Deep Ecliptic Survey team to compute a debiased 
//-magnitude distribution (Adams et ah, 2014). Fortunately, Centaurs have 


short dynamical lifetimes, so their size distribution is a relic of their source 
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population. The cold Kuiper Belt (e.g. Holman and Wisdom, 1993 Levison 


and Duncan,1997 1 

%aser et al. 

2010 Volk and Malhotra 

2011), the Neptune 

Trojans ( 

dorner and Lykawka 

2010), the inner Oort cloud (Emel’yanenko et 

al., 20051 

Kaib et al 

.1 20091 

Brasser et al., 

20121 Volk and Malhotra, 

20131 de 

la Fuente Marcos and de la Fuente Marcos 

, 2014 

Fouchard et ah, 2014), the 

Plutinos 

Morbidell 

i, 1997 

di Sisto et al. 

2010) 

and the scattered disk (di 


Sisto and Brunini 

2007 

Volk and Malhotra 

2008 

2013 

) could all be Centaur 

sources. However, 

no empirical information exists on the size distribution of 


objects in the Oort cloud, and Fraser et al. (2010) find that the scattered disk 


is not populous enough to explain the observed infiux of comets into the inner 


Solar System. Doressoundiram et al. (2005) also show that Centaur colors 


are not consistent with an origin in the scattered disk. Moreover, Schlichting 


et al. (2013) show that the cold Kuiper Belt and scattered disk objects have 


size distributions that follow the same functional form, only with different 
maximum sizes. Calculations presented here are based on the cold Kuiper 


Belt size spectrum of Schlichting et al. (2013), which is a close match to 


the size spectrum of Saturnian satellite impactors inferred from the cratering 


record (Minton et ah, 2012). Schlichting et al. (2013) used a combination 


of theoretical coagulation models, occultation surveys, and observations of 
large KBOs to constrain the size spectrum. 

The first estimate of the total number of Centaurs comes from the sim¬ 


ulations of Tiscareno and Malhotra (2003), who investigated the dynamical 


evolution of the observed Centaurs over 100 Myr. The top panel of Figure [T] 
shows their computed time-averaged eccentricity distribution. Tiscareno and| 
Malhotra (2003) also estimated the detection fraction of Centaurs as a func¬ 


tion of eccentricity, which is reproduced in the bottom panel of Figure The 
detection fraction estimate holds for Centaurs with i? > 30 km. Multiplying 
the eccentricity distribution with a fit to the detectability function (black 
line in the bottom panel of Figure and summing over the 0-1 eccentricity 
range yields an estimate of fdet = 3.7% for the fraction of Centaurs with 
/? > 30 km that have been detected. The total number of large Centaurs 
with i? > 30 km is then is ~ Nohs/fdet, where Whs = 53 is the number of 


Centaurs that had been discovered when the Tiscareno and Malhotra (2003) 
calculations were performed. 

The next step in determining the total number of Centaurs is to find the 
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radius of the largest Centaur. For R> 30 km, 


N>{R) = 


Nn 


c 


Ro 


i-C 


( 2 ) 


In E quation [2| A^>(i^inax) = 1, iV>(30 km) = A^obs//det = 1432 and C = 4 


(e.g. " 

?rujillo et ah, 2001 

Fraser et ah, 2008 [Minton et ah 

2012; 

Schlichting 

et ah 

2013 

), so that i?max = 338 km. An estimate of the total number 


of Centaurs then follows, given an analytical form for the differential size 


distribution dN/dR. Schlichting et ah (2013) hnd a KBO size distribution 
of the form dN/dR oc where = 2 for 10 km< R < 30 km; = 5.8 
for 2 km< i? < 10 km; and = 2.5 for 0.1 km< R < 2 km. We set a lower 
limit of i? = 1 km to the size of Centaurs considered here, which is justihed 
because the mass contained in the smallest bodies is negligible unless C ^ 4. 


The Schlichting et ah (2013) conclusion that C < 4 for the smallest bodies 


is supported by sky brightness measurements, which rule out C ^ 3.4 for 


7? < 1 km ( 

Kenyon and Windhorst 

2001; 

Ichikawa and Fukugita 

2011) 

The size distribution computed based on the 

Tiscareno and Malhotra 

(2003 


maximum Centaur-size estimate contains 2.8 x 10^ comets with R > 1 km, 
and is shown in the top pa nel of Figure (red curve). The distribution 
agrees well with the results of Sheppard et ah (2000), who predict about 100 


Centaurs with radii above 50 km. However, the number of small bodies is 


an order of magnitude lower than the di Sisto and Brunini (2007) estimate 
of ~ 2.8 X 10® Centaurs with radii above 1 km. 

Other estimates of the total number of Centaurs come from radius mea¬ 
surements of Centaurs and KBOs. The most conservative estimates come 
from assuming that Chariklo, the largest observed Centaur, is in fact the 
largest Centaur in the Solar System. (It is highly likely that the largest 
Centaur has not been observed, given that the detection probability is ex¬ 
tremely low for even moderately eccentric orbits.) Chariklo radius estimates 


range between 118 km and 151 km ( 

Fornasier et al. 

2013 

Stansberry et al. 

2008 

Groussin et al. 

2004; 

Altenhoff et al., 

2001 Jewitt and Kalas 

1998) 


The green and black curves in Figure show size distributions where the 
largest body takes on the maximum and minimum observational estimates 
of Chariklo’s radius, respectively. Finally, Figure shows a size distribution 
that is optimistic about the size of the largest body, with i?max = 458.5 km, 
the maximum measured radius of the Plutino Orcus (blue curve). Larger 
TNOs such as Quaoar and Pluto have higher densities that suggest differ- 
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entiation, whereas Orcus’ bulk density is more consistent with the undif¬ 
ferentiated comet population. It is plausible, then, that Orcus represents a 
transition object between Kuiper Belt comets/Centaurs and true dwarf plan¬ 
ets. Note that this “optimistic” size distribution predicts only a factor of 2.5 
more Centaurs than using the Tiscareno and Malhotra (2003) results, but 


brings the number of small bodies closer to the predictions of |di Sisto and 
Brunini ( 2007[ ). 


The bottom panel of Figure|^shows the cumulative mass function M<(i?), 
the mass of Centaurs with radii less than a given value. For each size distribu¬ 
tion from the top of Figure]^ three possible Centaur densities are considered. 
The highest density is the maximum inferred value for Orcus (1.53 g cm“^. 



Centaurs are given the same density regardless of size, though in reality the 
larger bodies probably have higher densities. 

It is important to note that the abundance of large Centaurs may be 
out of equilibrium with the source KBO population, due both to stochastic 
effects and to the possibility that the diffusion rate from the Kuiper Belt 
to the inner Solar System is a function of size. Since impacts by 100-km 
objects onto ice giants occur less than once per 10 Myr, these objects may be 
safely neglected from any discussion of the steady-state ice-giant luminosities. 
However, the size distributions presented here are anchored by the large 
objects, so disequilibrium from the source KBO population could lead to 
large uncertainties the number of small bodies. 


4. Impact-generated luminosity 

Given the Centaur size distribution, impact speed, and radiative time 
constants of the ice giant atmospheres, we can compute the impact-generated 
luminosity of ice giants as a function of time. The impact speed is = 
'^esc + (euorb)^, where Uesc is the planet’s escape speed, e is the Centaur’s 
eccentricity, and Uorb is the Centaur’s orbital speed. For Centaurs crossing 
the orbits of Uranus and Neptune, Uorb ^ and impacts occur roughly at 
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escape speed. Each impact converts Centaur kinetic energy into heat, which 
the ice giant atmospheres radiate away over a characteristic timescale. 

Radiative time constants Trad used here are 95 years for Uranus and 105 
years for Neptune. Conrath et ah (1990) computed radiative time constants 
based on prohles of Qir, the infrared energy deposited per unit volume and 
time at each level of the planet atmosphere. Qir is usually negative, indi¬ 
cating cooling, and is dominated by methane emission. To compute Trad, 


Conrath et al. (1990) hrst constructed model atmospheres by imposing the 


latitudnally dependent temperature measured by Voyager data as a bound¬ 
ary condition, then perturbed the temperature prohle at each atmosphere 
depth to compute the change in Qi^. Changes in the composition of ice 
giant atmospheres, such as aerosol darkening predicted to correspond with 


the 11-year solar cycle (Baines and Smith, 1990 Hammel and Lockwood 


2007), could alter the radiative time constants; how much is not clear as the 


Conrath et al. (1990) calculation has not been updated based on new data. 


Since radiative time constants become independent of atmospheric pressure 
for P > 0.5 bar (the boundary between the radiative and convective layers 
according to Conrath et al. (1990) and de Pater and Lissauer ( 2010[ )) and 
even the smallest impactors considered (i? = 1 km) would penetrate deeper 
than the 1-bar level, impact-generated luminosities are assumed not to be 
sensitive to the penetration depth of the impactors. 

We use a Monte Carlo simulation of successive ice giant impacts from 
Centaurs with radii chosen randomly from the size distributions plotted in 
Figure The Monte Carlo simulation reveals how often Centaur impacts can 
bring Neptune’s internal luminosity (excluding reprocessed energy from the 
Sun) up to the observed value of 3.64 x 10^^ erg s“^. A sequence of impacts 
meets the criteria for a Poisson process: gravitational interactions between 
Centaurs are negligible, so the probability of an impact is independent of 
previous impact history, and the likelihood of an impact in a given amount 
of time increases with the length of the time interval under consideration. 
Impacts are therefore modeled as a Poisson process with rate parameter 
A = l/xiinp, where Timp is the mean time interval between impacts, and the 
distribution of time intervals between impacts is exponential. Only impactors 
larger than 1 km are included in the simulation. Based on their models of 
diffusion from the Kuiper Belt, Levison and Duncan (1997) calculated that 
there should be 1.2 x 10^ ecliptic comets (Centaurs, Halley-type comets and 
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Jupiter-family comets combinecQ, assuming a maximum radius of -Rmax = 
200 km and a minimum radius -Rmin in the range 0.5-2 km. Levison and 


Duncan (1997) find impact rates of Aq = 1.2 x 10 ^ yr ^ for Neptune and 
Aq = 1.3 X 10“^ yr“^ for Uranus. Rate parameters are scaled to the Levison 


and Duncan (|1997|) results according to the number of Centaurs with radii 

iV>(l km) 


over 1 km, such that 


A — An 


1.2 X 10^ 


( 3 ) 


The simulation is built on the simplifying assumption that impact energy 
is deposited instantaneously and then redistributed throughout the atmo¬ 
sphere immediately. Subsequent to an impact, Newtonian cooling governs 
impact-generated planet luminosity as a function of time. Luminosity is then 
determined by both the impact energy Ui^p and the atmospheric radiative 
time constant: 


L = (Limp + Lo)e 

Limp L/imp/TVad 


( 4 ) 

( 5 ) 


In Equation timp is the time of the most recent impact and Lq represents 
any residual luminosity left over from previous impacts. We assume 100% 
conversion of impactor kinetic energy to heat, though in reality some energy 
will be lost to sublimation and ablation of the impactor. As the goal of 
this experiment is to assess how much of the ice giant energy balance can be 
explained by impacts, no other source of internal energy (such as gravitational 
contraction or radioactive decay) is included. 

The simulation proceeds as follows: (1) a random set of impact time in¬ 
tervals is drawn from an exponential distribution. (2) For each impact, a 
Centaur radius is randomly selected from the size distribution under consid¬ 
eration. Figure shows two realizations of randomly chosen impactor size 


^The size distributions based on measurements of Chariklo and Orcus do not encode 
any dynamical information about the impactors, so that they may represent Jupiter-family 
comets or Halley-type comets as well as Centaurs. No correction needs to be made to the 
total number of objects to be consistent with the Levison and Duncan (19971 impact 
rates. Of the estimates of the total number of Uranus and Neptune impactors in Section 
1^ only the one based on the Tiscareno and Malhotra (20031 detectability simulation is 
restricted to true Centaurs with Q > 30.2 AU. However, since Levison and Duncan ( |1997 ) 
find that true Centaurs are far more common than shorter-period comets, we also make 
no correction to the total number of impactors from this size distribution. 




























distributions. (3) The time evolution of the planet’s luminosity is calcu¬ 
lated. At each impact time, the luminosity Limp is added to any remaining 
luminosity from previous impacts Lq and a Newtonian cooling process starts 
(Equations 1^ and 1^. (4) The cumulative probability distribution of planet 
luminosity is calculated. 


5. Results and Conclusions 


Figure shows probability distributions of planet luminosity from a 10- 
Myr simulation, for all four size distributions shown in Figure and for 
Centaur densities of 0.5 g cm“^ (left) and 1.5 g cm“^ (right). For only the 
two most populous size distributions did the simulations record impacts that 
produce luminosities of order Neptune’s observed value, 3.64 x 10^^ erg s“^— 
and then extremely infrequently. If the total Centaur mass is of order Pluto’s 
mass (see blue dashed curve in the bottom panel of Figure [^, an impact that 
yields a luminosity of Ljmp > 10^^ erg s“^ occurs every 0.6 Myr. The largest 
luminosity generated by a single impact that the simulations produced was 
Limp = 3.3 X 10^^ erg s“^, from a 13-km impactor with p = 1.5 g cm“^. 

Based on current understanding of Centaur dynamics and radii, it is clear 
that impacts cannot explain the ice giants’ energy balance. Not only are 
energetic impacts far too rare to provide Neptune’s luminosity, but [Levison 
and Duncan (1997) show that Uranus and Neptune experience impacts at 


almost the same rate. Furthermore, their atmospheres radiate away energy 
on similar timescales. Only after a rare giant impact onto one planet are their 
impact-induced luminosities are substantially different. Instead, Uranus and 
Neptune’s different energy balances are probably the result of a structural 


dichotomy between the two planets (Nettelmann et ah, 2013), though shape 


and rotation data are not yet precise enough to elucidate exactly how the ice 
giant interior structures differ. 

There are two circumstances in which it may be worth revisiting the con¬ 
nection between impacts and ice giant energy balance. First, a re-calculation 
is warranted if observers discover a signihcant population of Neptune-crossing 
comets that have not been well represented in dynamical simulations (i.e. 
high-inclination objects). With a factor of 100 increase in Centaur abun¬ 
dance over and above the most populous Centaur size distribution consid¬ 
ered here, impacts would contribute substantially to Neptune’s steady-state 
energy balance. Furthermore, high-inclination Centaurs from the Oort cloud 


have longer dynamical lifetimes than those with low inclinations (e.g. di Sisto 
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and Brunini, 2007 Brasser et al., 2012; Volk and Malhotra, 2013), which may 


indicate that this work nnderestimates the steady-state Centanr abnndance. 
The nnmber of low-inclination Centanrs sonrced from the scattered disk and 
cold Kniper Belt is not likely to increase mnch as both popnlations have 


been well characterized (e.g. 

Adams et al., 

2014 

)• 

The highest estimate of 

the number of Centaurs (di Sisto and Brunini 

2007 

) is a mere factor of four 


higher than the maximnm nnmber considered here, and Horner et al. 
predict only 4.4 x 10^ Centanrs larger than 


( |2004[ ) 

1 km. We do not consider a 
snbstantial npward revision to the total nnmber of Centanrs likely. Second, 
if models of the KBO size distribntion are revised to inclnde a signihcantly 
higher proportion of ~ 10-km objects, median impact energies wonld pro- 
dnce Inminosities of order Neptnne’s cnrrent valne. However, for impact rates 
calcnlated by Levison and Dnncan (1997), changing the KBO size spectrnm 
to favor 10-km objects conld pnsh Uranns’ Inminosity well over observed 
valnes. In any case, planetesimal-growth simnlations, occnltation statistics, 
and crater size measnrements are converging in a way that make snbstantial 


revisions to the 10-km KBO abundance unlikely (e.g. 

Schlichting et al. 

2009 

Kenyon and Bromley 

2012; 

Minton et al. 

2012| Schlichting et al. 

2013 

)• 


One possibly frnitfnl avenne for fntnre research is investigating the con¬ 
nection between impacts and seismic distnrbances or storms. A comet of 
radius 1 km would deposit 1.2 x 10^® erg of energy locally over a timescale 
of ~ 1000 seconds, based on the Shoemaker-Levy Jupiter impact models of 


Harrington and Doming (2001). Marley (1994) demonstrated that impacts 


of energies 10^* ergs and higher would produce detectable temperature fluc¬ 
tuations exceeding 1 K on Jupiter due to the excitation of seismic waves. 
However, impacts by 1 km objects are still too infrequent—once per hundred 
years, even for generous estimates of the number of Centaurs such as those 
of di Sisto and Brunini (2007)—to make ice-giant seismology a reasonable 
observing campaign. Detection of seismic activity would be serendipitous. 
It is faintly possible, however, that sedimentation of particles ablated from 
an impactor could be observable. In Keck NIRC2 adaptive optics observa¬ 
tions of Uranus on 12 August 2004, Sromosky and Fry (2005) noted three 
low surface-brightness features that appeared and disappeared on timescales 


of 0.3-1.5 hours. Based on observations by Carlson et al. (1988), Sromosky 


and Fry (2005) suggest that sedimentation of 300-/im particles caused by 


vertical motion of a cloud could produce the rapid disappearance of the faint 
features. While vertical shear acting on clouds is the most likely reason for 
small particles to descend in Uranus’ atmosphere, we suggest cautiously con- 
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sidering impactors a source of settling particles. Extending the Schlichting et 


al. (2013) size distribution down to its lowest radius limit of 0.01 km yields 


a total of > 10^^ Centaurs, even assuming the smallest maximum Centaur 
size considered here (118 km; minimum measured Chariklo radius). At min¬ 
imum, tens of 0.01 km-size objects would impact each ice giant per year. 


Given that Uranus has recently shown record-breaking storms (de Pater et 


al., 2015), despite its equinox passage possibly decreasing convective activity. 


it is worth considering whether impacts may trigger local storms. 
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Centaur Detectability 
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Figure 1: Top: Time-average of Centaurs per eccentricity bin from the 100 Myr Centaur 
dynamics simulation of |Tiscareno and Malhotra (2003) (a reproduction of their Fig. 10a). 
Bottom: Centaur detection probability based on the Tiscareno and Malhotra ( |2003 ) 
Monte Carlo simulation (reproduction of their Fig. 10b). The black line shows an expo¬ 
nential ht to their data. Combining the two curves suggests that < 4% of Centaurs have 
been detected. 
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Figure 2: Top: Size distribution N^{R), the number of Centaurs with radii greater than 
size R. Colors show size distributions corresponding to different assumed sizes of the 
largest Centaur. For the red curve marked TM03, the largest Centaur radius is calculated 
self-consistently using the Centaur detection probability of Tiscareno and Malhotra (2003), 
which suggests that 3.7% of the Centaurs with i? > 30 km have been detected. The blue 
curve marked “Orcus” shows a size distribution for which the largest Centaur has the 
diameter of Orcus, 917 km. The green and black lines labeled Chariklo 1 and Chariklo 2 
show size distributions with the largest Centaur diameter given by the minimum and 
maximum measured values for Chariklo. Bottom: Mass distribution the total 

mass of Centaurs with radii less than size R. The color scheme is as above. Line styles 
denote assumptions about Centaur density. Orcus has a maximum measured density of 
1.5 g cm“^, Salacia has a maximum measured density of 1.3 g cm“^, and comet 67P 
(Churyumov-Gerasimenko) has a density of 0.5 g cm“^. 
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Figure 3: Sample size distributions of impactors chosen from 1 Myr Monte Carlo simula¬ 
tion. Stochastic behavior appears for larger Centaur radii. 




Figure 4: Cumulative probability distribution of impact-generated planet luminosity for 
each of the four size distributions under consideration. In the left panel, the model Cen¬ 
taurs have the same density as Comet 67P, p = 0.5 g cm“^. In the right panel, the 
assumed Centaur density is p = 1.5 g cm“^. The color scheme is the same as in Figure 
Solid lines show Neptune luminosities and dashed lines show Uranus luminosities. The 
minimum luminosity produced by an impact is Limp = 2 x 10^® erg s“^. Lower values 
indicates long time intervals between impacts. 
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